function [a,da,dda] = func_InteNewmark(c0,c1,c2,c3,c4,c5,c6,c7,a,da,dda,Q,M,C,KK)
%myFun - Description
%
% Syntax: [a,da,dda] = func_InteNewmark(a,da,dda,M,K,C,c0,c1,c2,c3,c4,c5,c6,c7)
%
% Long description

    QQ=Q+M*(c0*a+c2*da+c3*dda)+C*(c1*a+c4*da+c5*dda); % 如果是时间变化的外力，此处需修改
    
    af=a; % the former a
    daf=da;
    ddaf=dda;
    
    a=KK\QQ;
    
    dda=c0*(a-af)-c2*daf-c3*ddaf;
    da=daf+c6*ddaf+c7*dda; 
end